Slow light with three-level atoms in metamaterial waveguides 
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Metamaterial is promising for enchancing the capability of plasmonic devices. We consider a 
cylindrical waveguide with three- level A atoms embedded in the dielectric core. By comparing 
metal cladding vs metamaterial cladding of a waveguide with A atoms in the core, we show that, for 
a fixed amount of slowing of light due to electromagnetically induced transparency, the metamaterial 
cladding outperforms in terms of the inherent loss. 



I. INTRODUCTION 

Plasmonic devices show promise in increasing the speed of electronic devices and networks while still meeting 
the size requirements of modern electronic devices due to their ability to confine light to sub- wavelength scales jl| . 
Combining such devices with optical phenomena, such as slow light, could yield an array of new optical devices 
to augment or replace existing technology. The feasibility of using nonlinear metamaterials for slow light has been 
considered 0, Q , and a scheme for photon-echo quantum memory for surface plasmon-polaritons on a metamaterial 
interface has been proposed f^. Controllable slow light is useful for optical delay lines, optical buffers [HI and 
enhanced nonlinear interactions 6]; however, achieving highly-confined slow light with weak fields is impractical due 
to the high losses suffered as a result of the plasmonic confinement of light 7] . There exists a scheme @, using a 
fiat interface metamaterial-dielectric waveguide for slowing highly-confined light through electromagnetically induced 
transparency (EIT) while using metamaterial to minimize losses. EIT with weak fields is interesting in its own 
right for fundamental reasons, to distinguish it from Autler-Townes splitting, and also for applications to weak-field 
sensing However, the scheme has the limitation that field confinement in the waveguide is only in one transverse 
spatial dimension, which allows propagating fields to diverge. 

We propose using a cylindrical waveguide structure as shown in Fig. [1] rather than a fiat interface, to inherently 
provide confinement in both transverse directions and prevent field divergence. The waveguide is composed of a 
dielectric core and a metamaterial cladding that surrounds the core. Three-level A atoms (Fig. [5]) are homogeneously 
embedded throughout the core of the waveguide and are driven by a pump field in a low- loss surface mode [9[ , thereby 
enabling EIT for slowing the signal field in another low-loss surface mode. To prepare EIT the atoms are initially 
strongly pumped to depopulate the |s) level and relax to the |g) state. This process creates Fano interference in 
the |e) O |g) transition resulting in a transparency window at the |g) O |e) transition frequency. As a result of 
the refractive index change, the group velocity of the signal pulse is reduced when traveling through the transparent 
medium. After preparing for the EIT, the intensity of the pump field is lowered because the weaker the pump field is 
the more reduction of group velocity of the signal field. Hence, coherently controlled slow light is achievable by tuning 
the intensity of the pump field. Furthermore, by comparing the metamaterial-cladded waveguide to a metal-cladded 
one with the same permittivity, we show that a low-loss surface mode [12] of the metamaterial-cladded guide enables 
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FIG. 1. (a) The cylindrical waveguide cross section, (b) A schematic of the metamaterial-dielectric waveguide embedded with 
three-level A atoms (black dots) in the core. The outer layer (cladding) is a metamaterial whereas the core is a dielectric. 
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FIG. 2. A schematic of the three-level A atom. The ground state is jg), the excited state is |e) and |s) is a metastable state. 
The frequency of the pump field is ujp and of the signal field is ljs, and A and Ar are detunings. 



the same degree of slowing of light but with reduced losses. Our analysis is important to fabricating coherently- 
controlled low-loss slow light devices with current metamaterial technology and we provide conditions necessary for 
the existence of the low-loss surface mode for slowing light in the waveguide. 

In sectionini we develop the theory for field propagation in the cylindrical metamaterial-cladded waveguides without 
the doping of the three-level atoms. Then we extend the theory to allow for doping the atoms in the core of the 
waveguide. In section IIIIl we describe the numerical method we use to solve Maxwell's equations in the doped 
waveguide and present results obtained from the solutions. A discussion of our scheme is given in section ITVl and a 
summary of the work is given in section |Vl 



II. THEORY 



Metamaterials are designed with artificial structure, giving them electromagnetic responses different to those of their 
constituent materials. We consider a fishnet-type metamaterial as they have been well studied, both experimentally 
and theoretically [l3l - [l6| . Metamaterials of this design are inherently lossy so the permittivity and permeability 
expressions must necessarily include dissipation terms. The permittivity of the fishnet metamaterial is described by 
the Drude model [l7j . 

eMM(w) 

^ '^b -, -^TTTi (1) 

Co Uj{uj + lle) 

with eg the permittivity of free space, eb the background permittivity, lOc the plasma frequency of the metal, and Fo 
the electric damping term. The permeability, which is in practice achieved through the structure of the metamaterial, 
is described by the relation [isj 

= Mb + — 7 — --TT^i l^j 

Mo UJ^ - UJ [UJ + iTi-a) 

with fj,Q the permeability of free space, Mb the background permeability, F a geometric parameter, ujq the resonance 
frequency, and Fm the magnetic damping term. Equations ([IJ and ^ describe the electromagnetic properties of the 
metamaterial cladding and are used, along with a constant permittivity for the core, to determine the modes of the 
undoped metamaterial-dielectric waveguide. 

After preparation for EIT, the |s) state is not populated and the atoms are no longer excited by the pump field, 
which means the pump field effectively travels in the modes of the undoped waveguide. On the other hand, the signal 
field can excite the atoms, so the signal field travels in the modes of the doped waveguide. In the following, we develop 
the theory for finding the modes in the undoped waveguide before finding the modes of the doped waveguide. 

The guided modes in a cylindrical waveguide are transverse electric (TE), transverse magnetic (TM) or a combi- 
nation of the two (denoted HE or EH) [i9[. The TM and TE modes are those with field components Hz = and 
Ez ~ Q respectively, with z the propagation direction. We restrict our analysis to TM modes as they are invariant in 
azimuthal angle 0, which allows us to simplify the calculations. In general, the analysis can be done with any choice 
of modes. 
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For the TM modes of the metamaterial-dielectric waveguide, the electric fields have the form 

E{r, 0, z, t) = E{r, 0) e'CP^-'^t) ^ 

with c.c. indicating the complex conjugate, and j3 — P + ia the complex propagation constants for the allowed modes 
of the waveguide, which are obtained by solving the dispersion relation. Using the definitions 



dr 
and 

dKm in) 



K'm («7) 



dr 



(4) 



(5) 



with Jm and Km the Bessel function and modified Bessel function respectively, the dispersion relation for TM modes 
in an undoped cylindrical metamaterial-dielectric waveguide is [l9| 

Ed JL (aK) , Emm K'm (aj) ^ ^ 



k2 (an) 7^ Km (07) 



with 



7 := - w^eMMMMM, K := sfuS^^^jM--^ , (7) 



Ed the permittivity of the dielectric core, and a the core radius. With this dispersion relation, we can find the low loss 
TM surface mode for the pump field as detailed in (T2j . 

The susceptibility of the undoped waveguide core is simply the dielectric susceptibility, 

Xd = Ed/eo - 1- (8) 

When the core is doped with the three-level atoms, the susceptibility of the core becomes the sum Xd + x('^: which 
alters the guided modes near the EIT resonance frequency from those of the undoped guide. To find the low loss TM 
surface mode of the signal field, we must solve the wave equation with the susceptibility of the doped waveguide core. 
The susceptibility of the pumped three-level A atoms under EIT is [lo| 

x(w,r) = ^^ — , ^r~r^' 

nd^ 7eg-zA+ |17p(r)| (7sg-iAR) 

with c the speed of light in vacuum, — c^/edfM the refractive index of the dielectric core with /id the permeability 
of the dielectric core, 7ij the decay rate from level i to level j and 

flo 2 ~ 2 n Pa (10) 

with Pa the number density of three-level atoms. The approximation in Eq. (jlOp is valid near the EIT resonance 
frequency. Defining ujij as the transition frequency for the \i) O \j) transition the detunings are A ~ Wcg — Wg, 
Ar := a — Wos + i^p- The pump is set on resonance with the |s) O |e) transition such that Ar = A. The Rabi 
frequency of the pump is 

np{r)^^d-E^{r), (11) 

with d the dipole moment of the |s) o |e) transition. 

Equations ([9|) and ()lip have r dependence but no or z dependence. This is because the pump is in a TM 
mode, which has a radial dependence but no azimuthal dependence, implying (j) invariance. Furthermore, the pump 
is approximated as non-depleting, which implies z invariance. This approximation is valid for short propagation 
lengths, such that the pump does not deplete considerably. The results presented in this paper are based on this 
approximation. 

For longer propagation lengths the pump does deplete, which reduces the Rabi frequency flp and leads to a further 
reduction in the group velocity. In this case, one can treat the whole waveguide as a series of concatenated short 
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sections along the propagation direction, such that the pump intensity in each section is approximately constant. This 
treatment allows one to determine the group velocity of the signal in each section and calculate the overall delay of 
the signal in the whole waveguide. 

Inside the waveguide core, the pump field expression is given by 

E^(r) = A ^Jo(Kpr) z - t^JiM (12) 

with z and r unit vectors in the propagation and radial directions respectively, A a constant that determines the 
amplitude of the pump, /3p the complex propagation constant for the pump mode, and 



«P -\/w^^dMd^-^- (13) 
We define rj := Ad, which is a parameter adjusted by the amplitude of the pump, with d = d ■ z — d ■ r such that, 

n^ir) = I (^JoM - i^Ji(Kpr)j . (14) 

The permittivity of the doped core depends on the pump field. The pump field is in a TM surface mode (see Fig. [3]) , 
so the intensity varies steeply in the radial direction; thus, the signal field experiences a steeply graded refractive index 
in the core. As a result, within the core, the scalar wave equation used to solve for the modes of the undoped guide is 
no longer valid for the modes of the doped guide. In the cladding region, however, the permittivity and permeability 
are not spatially dependent, so we can still find the fields in this region using the scalar wave equation. 

To solve for the propagation constants and the fields in the doped guide, we must begin with the vector wave 
equation 

V^E + V ( ^-^E ■ Vecff(w, r) ) + k^Ceftiuj, r)E = 0, (15) 



with the Laplacian and 



e^^{u,r)^l ^o(l + Xd + x(-,r)), r < a. 



the effective permittivity of the waveguide with three-level A-atoms embedded in the core. As we are restricting our 
analysis to TM modes for both the pump and signal fields, Eq. can be decoupled into scalar equations with the 
z component taking the form 

T^^- + I - - 2 , f , eeff(w,r-) ) ^E, + Kts{uj,r)E,^ = 0, (17) 

dr^ \ r K^ff(cj,r)ecff(w,r) dr 
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with 

Keff(^, r) I \/^'Moecff(^,0-^', ^ < a, (18) 
\ ill r > a. 

Equation (jl7p is not generally analytically solvable, except at certain frequencies, so we numerically determine the 
dispersion and attenuation of the low-loss surface mode near the EIT resonance frequency. 



III. METHOD AND RESULTS 



There are two main problems with solving E g. (II 7p . First, j3 is complex, meaning it has two free parameters so a 
standard shooting method cannot be employed [20|. The second problem is that the fitness function. 



/ 



^0,clad 



(19) 



with H^ cmc and H^ ^iad the (j) component of the magnetic field in the core and cladding regions respectively, is not 
convex, so implementing a hill descent algorithm is not sufficient to determine if a particular value of /3 is a solution. 
To circumvent these problems we draw inspiration from the shooting method and employ a trial method to solve 
Eq. ^ for E,Xr,oj). 

We need to find Ez{r,uj) for a range of frequecies of interest. We begin at a frequency where an exact solution for 
a TM surface mode is found. To find the solution at an adjacent frequency, we slightly deviate the known /3 at the 
former frequency to obtain a trial value for /? and numerically solve Eq. (|17p for Ez{r,uj)] other components of the 
field are calculated from this trial Ez{r,uj). To test whether the trial /? is acceptable, we check if the corresponding 
trial Ez{r,ui) satisfies the boundary condition at the core-cladding interface to within a tolerance. Mathematically, 
the acceptance criterion is / < 10~^. If the trial -^^(r, w) does not satisfy the boundary condition, a different trial [3 
is chosen by slightly deviating the known f3 at the former frequency and the search is repeated. When an acceptable 
/? is found, we repeat the search for the next frequency until Ez{r,Lu) of the TM surface mode is determined for a 
range of frequecies of interest. 

We apply the same metamaterial parameters presented in our previous work (l2j to our calculations in this work. 
Hence, the TM surface mode of the undoped guide found in |12l | is the mode for the pump of our scheme. The 
metamaterial parameters are shown in Table HI The core parameters are Cd — 1.3eo, /^d — /^o a-nd the core radius 
a = Attc/uJc- For the three-level atoms embedded in the core, we set the number density = 1.26 x lO^^m"^ so 
that the resulting permittivity remains positive. The decay rate of the excited state is 7cg = 10^s~^ and that of the 
hyperfine state is 7sg = lO^^s^^; both of which are for Pr^^ 21 1. The mean signal frequency is Wg = i^og = 0.409a;e 



and the mean pump frequency is ujp = — 0.41cjc- The group velocity of the pump mode Vp = 0.47c. 



TABLE I. Parameter values for the metamaterial cladding. 



Parameter 


Value 


fb 


1 


Mb 


1 


CJc 


1.37 X lO^'s"' 


Fe 


2.73 X 10"s~' 


F^ 


Fc 


LUa 


0.2wc 


F 


0.5 



Using the method and parameters outlined above, we have calculated the propagation constants of the signal and 
pump TM surface modes of our scheme. With the complex propagation constant f3 = (3 + ia for a range of frequencies, 
we have calculated the group velocity of the signal field using Vg = (d/3/da;)~^. Figure S] shows the group velocity 
of the signal field surface mode in both the metal- and metamaterial-cladded waveguides for a range of rj. The plot 
shows that the group velocit;^f the signal for a metamaterial-cladded guide is proportional to 77^, and as rj oc |ri(r)|, 
this agrees with EIT theory TO]. 

Using the imaginary part of the propagation constant, we have obtained the attenuation curves for both the 
metamaterial- and metal-cladded waveguides as a function of detuning A as shown in Fig. [5] At the EIT resonance 
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FIG. 4. A plot of the numerically calculated group velocity at A = as a function of pump intensity for a metamaterial-cladded 
guide (O) with Fm ~ Fe/S and a metal-cladded guide (x). The solid curve is a least squares quadratic fit. 
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FIG. 5. Plots of the relative attenuation of the surface mode of the metamaterial guide to that of a metal guide (QMM/aMctai) 
with 7)/ft = 5 X 10*s"\ for Fm = Fc/3 (solid) and Fm = Fc/lOO (dashed). 




(i.e. zero detuning), the metamaterial-dielectric guide provides nearly a 20% reduction in attenuation over a metal- 
dielectric guide for Fm = Fo/S (solid line). For Fm = Fc/lOO (dashed line), the attenuation reduction improves 
to near 40% over a metal-dielectric guide. When F,„ = Fo the attenuation of a metamaterial-dielectric guide is 
comparable to a metal-dielectric guide. For layered structures, such as the fishnet design, previous work [l8| has 
shown that Fm < Fc and a value of Fm = Fo/S should be feasible with current metamaterial technology [131 ■ Hence, 
our cylindrical metamaterial-dielectric waveguide is capable of supporting a low loss TM surface mode. 

To understand why the attenuation of the low- loss surface mode does not approach zero near Wg, as it does for 
the flat interface we plot the fraction of the total energy of the mode that resides in both the core and cladding. 
Figure [S] shows the fractional energy, ^, of the TM surface mode of the cylindrical metamaterial-dielectric guide as a 
function of frequency. For a given mode the fractional energy is the ratio of the energy in one part of the guide, the 
core for instance, to the total energy of the mode. The solid line is for the fractional energy in the core whereas the 
dashed line is for the cladding. The fraction of energy in the dielectric does not reach one and some of the signal field 
penetrates into the metamaterial, leading to a non-zero attenuation. 
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FIG. 6. Plots of the fraction of the total energy in the core (solid) and the cladding (dashed) as a function of frequency for the 
low-loss surface mode of the cylindrical metamaterial-dielectric waveguide. 

The results above show that embedding three-level A atoms in the core of a metamaterial-dielectric waveguide 
provides coherently controllable slowing of signal field with reduced losses. In the next section, we discuss the loss 
mechanism in the metamaterial waveguide and explain how the existence of the low loss TM surface mode is due to 
Fm taking a value less than Fe . 



In our calculation, we have treated the claddings as if they are infinitely thick. In practice, the TM surface modes 
for the pump and signal fields only penetrate into the metamaterial to some finite depth, so a finite thickness of the 
cladding regions is sufficient. The skin depth for both the signal and pump fields is about 0.1 As with As the wavelength 
of the signal at the EIT resonance. As typical fishnet metamaterial has thickness of O.lAg [15], a few layers of fishnet 
metamaterial (with dielectric spacers in between) is sufficient for the cladding. 

We find that the attenuation dip near (Fig. [5]) is not as pronounced in the presence of transverse field con- 
finement in the cylindrical waveguide. The flat metamaterial-dielectric interface has near-zero attenuation near this 
frequency, whereas the attenuation dip in the cylindrical guide does not approach zero. The cause for the non-zero 
attenuation is that the transverse confinement limits the amount of electromagnetic energy that is ejected from the 
lossy metamaterial. 

To illustrate why transverse confinement lessens the attenuation dip near Wg, consider a cylindrical guide in the 
limit of an infinite radius. In this limit, the interface that the surface wave propagates along becomes fiat and the 
dispersion relation for the cylindrical waveguide reduces to that for a flat interface waveguide. In the case of the flat 
interface, the energy can reside almost completely on the dielectric side leading to near-zero loss @. Alternatively, 
for a cylindrical guide in the zero core radius limit, the core almost vanishes; hence, the energy propagates almost 
entirely in the lossy metamaterial, which leads to significant loss. For a waveguide with a finite radius, the fraction 
of the total energy that can be confined to the dielectric lies between and 1; thus, there is a non- negligible amount 
of attenuation due to the field interacting with the metamaterial cladding. 

We have considered fishnet metamaterial for our scheme but the scheme is not restricted to using fishnet meta- 
material. Other types of metamaterial are applicable, as long as they have permittivity and permeability functions 
resembling equations ^ and ([2|), with e'y^yi < 0, < < 1 and Fm < Fc. Using other metamaterials is of interest 
because when layered fishnet metamaterials are bent to form a cylindrical cladding, the metamaterial may not interact 
with the fields isotropically. There now exist various promising fabrication techniques to building bulk metamaterials 



We remark that our scheme is limited to narrowband signal pulses. This is because the spectral width of the EIT 
transparency window is related to the intensity of the pump field [lo| . As slowing light requires a low intensity pump 
field, the resultant transparency window is spectrally narrow. The signal pulse must fit within the transparent spectral 



IV. DISCUSSION 
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window, so the signal pulse must be narrowband. 



V. SUMMARY 



We have considered two cylindrical dielectric-core waveguides of equal size and dimensions, with one having a metal 
cladding and the other a fishnet metamaterial cladding; other types of metamaterial cladding are also considerable. 
The permittivities of both claddings are given by the Drude model, but the metal cladding has constant permeability 
and the metamaterial cladding is given by the modified-Drude model. The signal field propagating through the 
waveguides are controUably slowed via electromagnetically induced transparency (EIT) due to the three-level A 
atoms homogeneously doped throughout the dielectric core of each waveguide. We have shown that the low-loss TM 
surface mode exists in the metamaterial waveguide provided that the permeability of the metamaterial cladding has 
a real part less than unity, i.e. Re (/u) < 1, and the magnetic damping rate is less than the electric damping rate, 
i. e. Fm < Fc. Previous work has shown that this condition is always true for certain optical metamaterials [18|. We 
predict that such cylindrical metamaterial-cladded waveguide is capable of delivering the same slowing of light but 
with reduced loss compared to the cylindrical metal-cladded waveguide. 
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